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Injecting spin-polarized carriers into semiconductor lasers provides important opportunities to 
extend what is known about spintronic devices, as well as to overcome many limitations of con¬ 
ventional (spin-unpolarized) lasers. By developing a microscopic model of spin-dependent optical 
gain derived from an accurate electronic structure in a quantum well-based laser, we study how 
its operation properties can be modified by spin-polarized carriers, carrier density, and resonant 
cavity design. We reveal that by applying a uniaxial strain, it is possible to attain a large birefrin¬ 
gence. While such birefringence is viewed as detrimental in conventional lasers, it could enable fast 
polarization oscillations of the emitted light in spin lasers which can be exploited for optical com¬ 
munication and high-performance interconnects. The resulting oscillation frequency (> 200 GHz) 
would significantly exceed the frequency range possible in conventional lasers. 

PACS numbers: 42.55.Px, 78.45. +h, 78.67.De, 78.67.Hc 


I. INTRODUCTION 

Both spin lasers and their conventional (spin- 
unpolarized) counterparts share three main elements: (i) 
the active (gain) region, responsible for optical amplifica¬ 
tion and stimulated emission, (ii) the resonant cavity, and 
(iii) the pump, which injects (optically or electrically) en- 
ergy/carriers. The main distinction of spin lasers is a net 
carrier spin polarization (spin imbalance) in the active 
region, which, in turn, can lead to crucial changes in 
their operation, as compared to their conventional coun¬ 
terparts. This spin imbalance is responsible for a cir¬ 
cularly polarized emitted light, a result of the conserva¬ 
tion of the total angular momentum during electron-hole 
recombination.— 

The experimental realization of spin lasers^^ii presents 
two important opportunities. The lasers provide a path 
to practical room-temperature spintronic devices with 
different operating principles, not limited to magnetore¬ 
sistive effects, which have enabled tremendous advances 
in magnetically stored informationi^^^l This requires re¬ 
visiting the common understanding of material parame¬ 
ters for desirable operation^ as well as a departure from 
more widely studied spintronic devices, where only one 
type of carrier (electrons) plays an active role. In con¬ 
trast, since semiconductor lasers are bipolar devices, a 
simultaneous description of electrons and holes is crucial. 

On the other hand, the interest in spin lasers is not 
limited to spintronics, as they may extend the limits of 
what is feasible with conventional semiconductor lasers. 
It was experimentally demonstrated that injecting spin- 
polarized carriers already leads to noticeable differences 
in the steady-state operation— The onset of lasing is 
attained for a smaller injection, lasing threshold reduc¬ 
tion, while the optical gain differs for different polariza¬ 
tions of light, leading to gain asymmetry, also referred to 
as gain anisotropy— In the stimulated emission, even 
a small carrier polarization in the active region can be 


greatly amplified and lead to the emission of completely 
circularly polarized emitted light, an example of a very 
efficient spin filtering—^ 



Figure 1: (Color online) Bucket model for (a) a conventional 
laser and (b) a spin laser.— Water added to the bucket rep¬ 
resents the carriers, and the water coming out represents the 
emitted light. Small leaks depict spontaneous emission, and 
overflowing water reaching the large opening corresponds to 
the lasing threshold. In (b) the two halves represent two spin 
populations (hot and cold water in the analogy), and they 
are filled separately. The partition between them is not per¬ 
fect: spin relaxation can cause the two populations to mix. 
The color code indicates conservation of angular momentum; 
an unpolarized pumping (violet) is an equal mixture of two 
polarized contributions (red and blue). 


An intuitive picture for a spin laser is provided by a 
bucket model in Fig.fTp^i^l The uneven water levels rep¬ 
resent the spin imbalance in the laser, which implies the 
following: (i) Lasing threshold reduction - in a parti¬ 
tioned bucket, less water needs to be pumped for it to 
overfill. There are also two thresholds (for cold and hot 
water) 1^ (ii) Gain asymmetry - an unequal amount of 
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hot and cold water comes out. A small spin imbalance of 
pumped carriers can (the two water levels slightly above 
and below the opening, respectively) result in a com¬ 
plete imbalance in the polarization of the emitted light 
(here only hot water gushes out) and, consequently, spin¬ 
filtering. These effects are attained at room temperature 
with either optical or electrical injection. The latter ex¬ 
perimental demonstratioiiii is a breakthrough towards 
practical use of spin lasers. 

Perhaps the most promising opportunity to overcome 
the limitations of conventional lasers lies in the dy¬ 
namic operation of spin lasers, predicted to provide en¬ 
hanced modulation bandwidth, improved switching prop¬ 
erties, and reduced parasitic frequency modulation, i. 
e., chirp Moreover, experiments have confirmed 
that in a given device a characteristic frequency of po¬ 
larization oscillations of the emitted light can signifi¬ 
cantly exceed the corresponding frequency of the inten¬ 
sity oscillationsThis behavior was attributed to 
birefringence - an anisotropy of the index of refraction, 
considered detrimental in conventional lasers^i 

What should we then require to attain high-frequency 
operation in spin lasers? Can we provide guidance for 
the design of an active region and a choice of the reso¬ 
nant cavity? Unfortunately, to address similar questions, 
we cannot simply rely on the widely used rate-equation 
description of spin lasers but instead we need 

to formulate a microscopic description. The crucial con¬ 
sideration is detailed knowledge of the spectral (energy- 
resolved) optical gain obtained from an accurate descrip¬ 
tion of the electronic structure in the active region, al¬ 
ready important to elucidate a steady-state operation of 
a spin laser. 

A typical vertical geometry, the so-called vertical cav¬ 
ity surface emitting lasers (VCSELs)^i^“— used in 
nearly all spin lasers, is illustrated in Fig. [Ha). Even 
among conventional lasers, VCSELs are recognized for 
their unique properties, making them particularly suit¬ 
able for optical data transmission)^ The resonant cavity 
is usually in the range of the emission wavelength, pro¬ 
viding a longitudinal single-mode operation. It is formed 
by a pair of parallel highly reflective mirrors made of 
distributed Bragg reflectors (DBRs), a layered structure 
with varying refractive index. The gain active (gain) re¬ 
gion, usually consists of III-V quantum wells (QWs) or 
quantum dots (QDs))^~— 

The key effect of the active region is producing a stim¬ 
ulated emission and coherent light that makes the laser 
such a unique light source. The corresponding opti¬ 
cal gain that describes stimulated emission, under suf¬ 
ficiently strong pumping/injection of carriers, can be il¬ 
lustrated pictorially in Figs.jHb) and IHc) for both con¬ 
ventional and spin lasers, respectively. In the latter case, 
it is convenient to decompose the photon density into 
different circular polarizations and distinguish that the 
gain is generally polarization-dependent. If we neglect 
any losses in the resonant cavity, such a gain would pro¬ 
vide an exponential growth rate with the distance across 
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Figure 2: (a) (Color online) Geometry of a vertical cavity 
surface-emitting laser. The resonant cavity of length L is 
formed between the two mirrors made of distributed Bragg 
reflectors (DBRs). The shaded region represents the active 
(gain) region of length Lw- The profile of a longitudinal op¬ 
tical mode is sketched. Schematic of the optical gain, g, in 
the active region for a conventional laser (b) and a spin laser 
(c). With external pumping/injection, a photon density S in¬ 
creases by 5S as it passes across the gain region.— In the spin 
laser, this increase depends on the positive (+)/negative( —) 
helicity of the light, S = S'*" + S~. 


a small segment of gain material.— Since both static and 
dynamic operations of spin lasers depend crucially on 
their corresponding optical gain, our focus will be to pro¬ 
vide its microscopic description derived from an accurate 
electronic structure of an active region. 

After this Introduction, in Sec. II we provide a theo¬ 
retical framework to calculate the gain in quantum well- 
based lasers. In Sec. Ill, we describe the correspond¬ 
ing electronic structure and the carrier populations un¬ 
der spin injection, the key prerequisites to understand¬ 
ing the spin-dependent gain and its spectral dependence, 
discussed in Sec. IV. Our gain calculations in Sec. V ex¬ 
plain how the steady-state properties of spin lasers can be 
modihed by spin-polarized carriers, carrier density, and 
resonant cavity design. In Sec. VI, we analyze the influ¬ 
ence of a uniaxial strain in the active region, which intro¬ 
duces a large birefringence with the resulting oscillation 
frequency that would significantly exceed the frequency 
range possible in conventional lasers. In Sec. VII, we de¬ 
scribe various considerations for the optimized design of 
spin lasers and the prospect of their ultrahigh-frequency 
operation. A brief summary in Sec. VIII ends our paper. 


II. THEORETICAL FRAMEWORK 


While both QWs and QDs^“— are used for the active 
region of spin lasers, we focus here on the QW implemen¬ 
tation also found in most of the commercial VCSELsi^ 
To obtain an accurate electronic structure in the active 
region, needed to calculate optical gain, we use the 8x8 
k-p method— The total Hamiltonian of the QW system, 
with the growth axis along the z direction, is 

Hqw{z) = Hkp{z) -b Hst{z) + Ho{z), 


( 1 ) 
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where Hkp{z) denotes the k-p term, Hst{z) describes the 
strain term, and Hq(z) includes the band-offset at the 
interface that generates the QW energy profile. The ex¬ 
plicit form of these different terms for zinc-blende crystals 
is given in Appendix A. 

Because common nonmagnetic semiconductors are well 
characterized by the vacuum permeability, po, a complex 
dielectric function e{ui) = er(uj) + £i(uj), where cu is the 
photon (angular) frequency, can be used to simply ex¬ 
press the dispersion and absorption of electromagnetic 
waves. The absorption coefficient describing gain or loss 
of the amplitude of an electromagnetic wave propagat¬ 
ing in such a medium is the negative value of the gain 
coefficient (or gain spectrum) 

5“(u;) = -—e“(a;), (2) 

CUr 


where c is the speed of light, is the dominant real part 
of the refractive index of the material;^ and e“(w) is the 
imaginary part of the dielectric function which generally 
depends on the polarization of light, a, given by 


I / 

c^v^k 


e“(w) = Co V - /J (5 - feu] , (3) 


where the indices c (not to be confused with the speed of 
light) and v label the conduction and valence subbands, 
respectively, k is the wave vector, - is the interband 
dipole transition amplitude, is the Fermi-Dirac dis¬ 
tribution for the electron occupancy in the conduction 
(valence) subbands, h is the Planck’s constant, is 
the interband transition frequency, and 6 is the Dirac 
delta-function, which is often replaced to include broad¬ 
ening effects for finite lifetimes.— The constant Co is 
Co = 47r^e^/(eomgtu^fl), where e is the electron charge, 
rriQ is the free electron mass, and D is the QW volume. 

Analogously to expressing the total photon density in 
Fig. m as the sum of different circular polarizations, 
s = s+ + S- ', in spin-resolved quantities we use sub¬ 
scripts to describe different spin projections, i. e., eigen¬ 
values of the (7z Pauli matrix. The total electron/hole 
density can be written as the sum of the spin up (-I-) and 
the spin down (—) electron/hole densities, n = n_|_ -|- n_ 
and p = p+ + P-. In this convention^i^i2£ using the 
conservation of angular momentum between carriers and 
photons, the recombination terms are n+p+, n_p_, while 
the corresponding polarization of the emitted light de¬ 
pends on the character of the valence band holes— We 
can simply define the carrier spin polarization 

Pa = {a+- a-)/{a++ a-), (4) 

where a = n, p— 

Using the dipole selection rules for the spin-conserving 
interband transitions, the gain spectrum. 


(5) 


can be expressed in terms of the contributions of spin up 
and down carriers. To obtain 5 “(_j(w), the summation 
of conduction and valence subbands is restricted to only 
one spin: ^ ^ Fr Eq. ([3|). 

C C-l-(-) V v + i-) 



Figure 3: (Color online) (a) Energy band diagram with a 
bandgap Eg and chemical potentials in conduction (valence) 
bands, p.c (pv) that in the presence of spin-polarized carri¬ 
ers become spin-dependent: pc{v)+ / Pciy)-- Unlike the 
rest of our analysis, here holes are spin-polarized, (b) Gain 
spectrum for unpolarized (solid) and spin-polarized electrons 
(dashed curves). Positive gain corresponds to the emission 
and negative gain to the absorption of photons. The gain 
threshold gth, required for lasing operation, is attained only 
for S~ helicity of light. 


To see how spin-polarized carriers could influence the 
gain, we show chemical potentials, pc{v)^ for a simplified 
conduction (valence) band in Fig. [31(a). The spin im¬ 
balance in the active region implies that pc(v) ''^Fl also 
become spin-dependent. Such different chemical poten¬ 
tials lead to the dependence of gain on the polarization of 
light, described in Fig.jSKb). Without spin-polarized car¬ 
riers, the gain is the same for S'"*" and S~ helicity of light. 
In an ideal semiconductor laser, g > 0 requires a popula¬ 
tion inversion for photon energies. Eg < Hu < {pc — k-v)- 
However, a gain broadening is inherent to lasers and, as 
depicted in Fig. IHb), g > 0 even below the bandgap, 
Hu < Eg. If we assume 0 [recall Eq. (I4|)] and 

Pp = 0 (accurately satisfied, as spins of holes relax much 
faster than electrons), we see different gain curves for 5+ 
and S~. The crossover from emission to absorption is 
now in the range of {pc- — gv-) and {pc+ — g-v+)- 
Optical injection of spin-polarized electrons is the most 
frequently used method to introduce spin-imbalance in 
lasers. Some spin lasers are therefore readily available 
since they can be based on commercial semiconductor 
lasers to which a source of circularly polarized light is 
added subsequently^ Such spin injection can be readily 
understood from dipole optical selection rules which ap¬ 
ply for both excitation and radiative recombination.— 

A simplified band diagram for a zinc-blende QW semi¬ 
conductor with the corresponding interband transitions 


3“(a;) = gl{u) + g'Liu) 
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is depicted in Fig.|4l At the Brillouin zone center, the va¬ 
lence band degeneracy of heavy and light holes (HH, LH) 
in the bulk semiconductor is lifted for QWs due to quan¬ 
tum confinement along the growth direction. The angu¬ 
lar momentum of absorbed circularly polarized light is 
transferred to the semiconductor. Electrons’ orbital mo¬ 
menta are directly oriented by light and, through spin- 
orbit interaction, their spins become polarizedii While 
initially holes are also polarized, their spin polarization is 
quickly lost^S Thus, as in Fig. [31[b), we assume through¬ 
out this work Pp = 0, unless stated otherwise. 
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Figure 4: (Color online) Schematic band structure and op¬ 
tical selection rules in zinc-blende QWs. (a) Conduction 
band (CB) and valence band with heavy and light holes (HH, 
LH) labeled by their total angular momentum j = 1/2 and 
j — 3/2, representing the states of the orbital angular momen¬ 
tum I = 0 and 1, respectively (Appendix A), (b) Interband 
dipole transitions near the band edge of a QW for light with 
positive and negative helicity, S^, between the sublevels la¬ 
beled by rrij, the projection of the total angular momentum 
on the -|- 2 :-axis (along the growth direction). The small ar¬ 
rows represent the projection of spin 1/2 of the orbital part 
that contributes to the transition, indicating that dipole tran¬ 
sitions do not change spin (Appendix B). 


Fig. [5l Our calculations, based on the k-p method and 
the 8x8 Hamiltonian from Eq. ([1]) (Appendix A), yield 
two confined CB subbands: CBl, CB2, and five VB sub¬ 
bands, labeled in Fig. [SKa) by the dominant component 
of the total envelope function at fc = 0. The larger num¬ 
ber of confined VB subbands stems from larger effective 
masses for holes than electrons4^ These differences in 
the effective masses also appear in the DOS shown in in 

Fig.EIb). 



Figure 5: (Color online) (a) Band structure for the 
Alo.sGao.rAs/GaAs QW for different fc-directions along [100], 
[010], and [110]. (b) DOS calculated from (a). The conduc¬ 
tion band DOS is multiplied by a factor of 20 to match the 
valence band scale. The bandgap is Eg — 1.479 eV (CBl-HHl 
energy difference). 


The spin polarization of excited electrons depends on 
the photon energy for S~^ or S~ light. From Fig.|31]b) we 
can infer that if only CB-HH are involved, |P„| —?► 1. At 
a larger Huj, involving also CB-LH transitions, |P„| is re¬ 
duced. Expressing oc where is the spherical 
harmonic, provides a simple connection between dipole 
selection rules and the conservation of angular momen¬ 
tum in optical transitions (Appendix B). 


III. ELECTRONIC STRUCTURE 

For our microscopic description of spin lasers we fo¬ 
cus, on a (Al,Ga)As/GaAs-based active region, a choice 
similar to many commercial VGSELs. We consider an 
Alo.sGao.rAs barrier and a single 8 nm thick GaAs QW4i 
The corresponding electronic structure of both the band 
dispersion and the density of states (DOS) is shown in 


As we seek to describe the gain spectrum in the ac¬ 
tive region, once we have the electronic structure, it is 
important to understand the effects associated with car¬ 
rier occupancies though injection/pumping [recall Fig. 2, 
Eqs. (2) and (3)]. In Figs. [6] (a), (c), and (e) we show 
both examples of injected unpolarized (P„ = 0) and spin- 
polarized {Pn = 0.5) electrons as seen in the equal and 
spin-split GB chemical potentials, respectively. The car¬ 
rier populations^ is given in Figs.[6l[b), (d), and (f) using 
the product of the Fermi-Dirac distribution and the DOS 
for CB and VB for both spin projections. 

IV. UNDERSTANDING THE 
SPIN-DEPENDENT GAIN 

From the conservation of angular momentum and 
polarization-dependent optical transitions, we can under¬ 
stand that even in conventional lasers, carrier spin plays 
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Figure 6: (Color online) Band structure of Fig. [SDa) with 
electron occupancy for (a) Pn=0, (c) P„ = 0.5, and (e) Pp — 
0. Carrier population expressed as a product of DOS from 
Fig-EIb) and the Fermi-Dirac distribution of electrons for (b) 
Prt=0, (d) Pn = 0.5, and (f) Pp — 0. The carrier density 
is fixed at n = p = 3 x 10^^ cm~^ and T — 300 K. The 
negative (positive) side of the x-axis represents spin down 
(up) electrons, dashed lines denote chemical potentials. The 
CB population is multiplied by 60 and shown in the same 
scale as for the VB. 


a role in determining the gain. However, in the absence 
of spin-polarized carriers^ the gain is identical for the 
two helicities: g'^ =5”, and we recover a simple descrip¬ 
tion (spin- and polarization-independent) from Fig. 2(b). 
In our notation, the upper (lower) index refers to the 
circular polarization (carrier spin) [recall Eq. (5)]. 

This behavior can be further understood from the 
gain spectrum in Figs. [7l[a) and (b), where we rec¬ 
ognize that g'^ = g~ requires: (i) g^ = g^ and 
g^ = gZ, dominated by CBl-HHl (1.479 eV = Eg) 
and CBl-LHl (1.501 eV) transitions, respectively (re¬ 
call Fig. [5]). No spin-imbalance implies spin-independent 
fic and gv [Fig. 3(a)] and thus and 

all vanish the photon energy Eph = huj = gc — 
fly- Throughout our calculations we choose a suitable 


Conventional laser Spin-laser 



Figure 7: (Color online) Gain spectra shown as a function of 
photon energy measured with respect to the energy bandgap. 
Conventional laser, = 0 for (a) and (b) S~ light po¬ 
larization. spin lasers, Pn ~ 0.5 for (c) and (d) S~ light 
polarizations. The carrier density n = p = 3 x 10^^ cm“^ and 
T — 300 K are the same as in Fig. (6] 


cosh~^ broadening^ with a full width at half-maximum 
(FWHM) of 19.75 meV, which accurately recovers the 
gain of conventional (Al,Ga)As/GaAs QW systems. 

We next turn to the gain spectrum of spin lasers. Why 
is their output different for S'"*" and S~ light, as depicted 
in Fig. 2(b)? Ghanging only = 0.5 from Figs.[7][a) and 
(b), we see very different results for 5'+ and S~ light in 
Figs.[7][c) and (d). > 0 implies that pLc+ > MC- [see 

Fig.IHc)], leading to a larger recombination between the 
spin up carriers (n+p+ > n-p-) and thus to a larger gy 
for 5"+ and S~ (red/dashed line) than g- (blue/dashed 
line). The combined effect of having spin-polarized car¬ 
riers and different polarization-dependent optical tran¬ 
sitions for spin up and down recombination is then re¬ 
sponsible for g'^ ^ g~, given by solid lines in Figs. [71[c) 
and[71[d). For this case, the emitted light S~ exceeds 
that with the opposite helicity. S'"*", and there is a gain 
asymmetryf^i^i^ another consequence of the polarization- 
dependent gain. The zero gain is attained at iJ-c+ — Pv 
for spin up (red curves) and pc- — Pv for spin down 
transitions (blue curves). The total gain, including both 
of these contributions, reaches zero at an intermediate 
value. Without any changes to the band structure, a sim¬ 
ple reversal of the carrier spin-polarization, > —Pm 
reverses the role of preferential light polarization. 
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V. STEADY-STATE GAIN PROPERTIES 

Within our framework, providing a spectral informa¬ 
tion for the gain, we can investigate how the carrier den¬ 
sity and its spin polarization, which can be readily mod¬ 
ified experimentally, can change the steady-state opera¬ 
tion of spin lasers. Specific to VCSELs, it is important to 
examine how their laser operation is related to the choice 
of a resonant cavity which defines the photon energy at 
which the constructive interference takes place. 



Eph - Eg (meV) 

Figure 8: (Color online) Evolution of the gain spectra with 
carrier density for: (a) = 0, (b) = 0.1, (c) Pn = 0.5, 

and (d) Pn = 1.0. In order to achieve emission, a certain 
value of carrier density should be added to the system. The 
second peak at Eph — Eg ~ 150 meV is related to transitions 
of CB2-HH2. Transitions of CB2-LH2 at Eph — Eg ~ 200 meV 
can be seen in the broader second peak for P„ = 1.0. The 
difference between g'^ and g~ that arise due to spin-polarized 
carriers in the system increases with Pn. For Pn = 1.0 there 
is no emission of 5"^ polarized light, i. e., this component is 
totally absorbed by the system. The diagonal arrow in Fig. 
|8] indicates the increase of carrier density in the curves. 


Most of the QW-based lasers do not have a doped 
active region, and they rely on a charge neutral car¬ 
rier injection (electrical or optical).— Here we choose 
n = p = 1,3, 5,7 X 10^^ cm“^, and spin polarizations 
Pn = 0,0.1,0.5,1, respectively. Electrical injection in in¬ 
trinsic III-V QWs using Fe or FeCo allows for |P„| ^ 
0.3 — 0.7^“— while |P„| —>■ 1 is attainable optically at 
room temperature — In most of the spin lasers, |P„| ^ 0.2 
in the active region. We focus on three resonant cavity 
positions: Ci, C 2 , C 3 (vertical lines), defining the corre¬ 
sponding energy of emitted photons Ci = 1.48 eV ~ 1.479 
eV (CBl-HHl transition), Ca = 1.5 eV ~ 1.501 eV (CBl- 
LHl transition), and C 3 = 1.52 eV (at the high energy 
side of the gain spectrum). 

The corresponding results are given in Fig. [ 8 ] showing 
gain spectra different for and S~. This gain asym¬ 
metry, ^ g~, is more pronounced at larger P„; at 
P„ = 1, there is even no emission. While this trend 
is expected and could be intuitively understood, there is 
a more complicated dependence of the gain asymmetry, 
g~ (tvjj) — g^{huj) on the carrier density and the choice 
of the detuning)^ the energy (frequency) difference be¬ 
tween the gain peak and the resonant cavity position. 



Figure 9: (Color online) Gain asymmetry obtained from Fig. [8] 
for: (a) P„ — 0.1, (b) Pn = 0.5, and (c) Pn = 1.0. As 
more carriers are added to the system, the asymmetry peak 
shifts to higher energies, however, this energy region is not 
necessarily in the regime of a positive gain. Gain asymmetry 
as a function of carrier density for: (d) Pn = 0.1, (e) Pn = 0.5, 
and (f) Pn = 1.0. Similar to the case of Figs. [IDa)-(c), the 
asymmetry peaks may not correspond to positive gain. 
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The gain asymmetry is one of the key figures of merit 
for spin lasers, and it can be viewed as crucial for their 
spin-selective properties, including robust spin-filtering 
or spin-amplification, in which even a small (few per¬ 
cent) in the active region leads to an almost complete 
polarization of the emitted light (of just one helicity)j^ 
Unfortunately, how to enhance the gain asymmetry, be¬ 
yond just increasing is largely unexplored. 



Figure 10: (Color online) Gain as a function of carrier density 
for: (a) P„ = 0.1, (b) P„ = 0.5, and (c) Pn — 1.0, with the 
cavity choices ci, C 2 , and cs. Comparing: (i) solid and short- 
dashed lines we can examine the spin-filtering effect, and (ii) 
solid and long-dashed curves, we can examine the threshold 
reduction. The solid horizontal line indicates the gain thresh¬ 
old, i.e., the losses in the cavity. To achieve the lasing, the 
value of gain must be greater than the gain threshold. 


To establish a more systematic understanding of a 
gain asymmetry, we closely examine g~{huj) — g'^(huj) in 
Figs.[ni(a)-inKc) for different Pn, carrier densities, and res¬ 
onant cavities. Increasing n, the anisotropy peak shifts 
to higher tvjj, indicating an occupation of higher energy 
subbands. However, the absolute anisotropy peak is not 
always in the emission region. For a desirable operation 
of a spin laser we should seek a large gain anisotropy 
with a positive (and a preferably large) gain. Comple¬ 


mentary information is given by Figs. |9Kd)-|9Kf) with a 
density evolution of g~ — g'^ for different cavity positions 
and Pn- Again, we see that the gain asymmetry peak 
can be attained outside of the lasing region. 

The results in Fig.|9]have shown a complex evolution of 
the gain asymmetry with the cavity position and carrier 
density. We now repeat a similar analysis for the gain 
itself in Fig. (TUI The gain calculated for two helicities and 
unpolarized light = S~), provides a useful guidance 
for the threshold reduction and the spin-filtering effect, 
invoked in a simple bucket model from Fig. [TJ 

We first consider = 0.1 which shows a behavior 
with an increase in n or, equivalently, an increase in in¬ 
jection, that could be expected from the bucket model. 
The threshold value of the gain (the onset of an over¬ 
flowing bucket), gth, is first reached for S~, then for un¬ 
polarized light, a sign of threshold reduction, and finally 
for S+ ( a subdominant helicity from the conservation of 
angular momentum and > 0). Therefore there is a 
spin-filtering interval of n (small, since itself is small) 
where we expect lasing with only one helicity. A similar 
behavior appears for all the cavity choices Ci, C 2 , and C 3 . 

We next turn to Pn = 0.5 where Ci shows trends ex¬ 
pected both from the bucket model and early work on 
spin lasersi^i^ An increase from = 0.1 to 0.5 enhances 
the threshold reduction and the spin-filtering interval. 
However, different cavity positions C 2 and C 3 reveal a dif¬ 
ferent behavior. There is a region where unpolarized light 
S~^ = S~ (long dashed lines) yields a greater gain than 
for S~ (solid lines). For C 3 the threshold is attained at 
smaller n for unpolarized light than for negative helicity, 
i.e., there is no threshold reductioiii^ With Pn = 1.0, 
the threshold reduction is only possible for Ci. 

These results reinforce the possibility for a versatile 
spin-VCSEL design by a careful choice of the resonant 
cavity, but they also caution that, depending on the given 
resonant cavity, the usual intuition about the influence 
of carrier density and spin polarization on the laser op¬ 
eration may not be appropriate. 


VI. STRAIN-INDUCED BIREFRINGENCE 

An important implication of an anisotropic dielectric 
function is the phenomenon of birefringence in which the 
refractive index, and thus the phase velocity of light, 
depends on the polarization of light Due to phase 
anisotropies in the laser cavity;^ the emitted frequen¬ 
cies of linearly polarized light in the x- and y-directions 
(S'® and S*') are usually different. Such birefringence is 
often undesired for the operation of conventional lasers 
since it is the origin for the typical complex polarization 
dynamics and chaotic polarization switching behavior in 
VCSELsj ^^:^^~ — While strong values of birefringence are 
usually considered to be an obstacle for the polarization 
control in spin-polarized lasers^^d^ the combination of a 
spin-induced gain asymmetry with birefringence in spin- 
VCSELs allows us to generate fast and controllable os- 












dilations between 5'+ and S~ polarizations j ^^ The fre¬ 
quency of these polarization oscillations is determined by 
the linear birefringence in the VCSEL cavity, and it can 
be much higher than the frequency of relaxation oscilla¬ 
tions of the carrier-photon system in conventional VC- 
SELs. This may pave the way towards ultrahigh band¬ 
width operation for optical communications 



k(nm''') k(nm''') 

Figure 11: (Color online) Band strncture with uniaxial strain 
in the active region for (a) ~ 0.019% and (b) e^x ~ 

0.058%. The inset shows a zoom around the HHl and LHl in¬ 
teraction region, where the difference between [100] and [010] 
directions is more visible. The energy gap of the system is 
Eg ~ 1.483 eV for case (a) and Eg ~ 1.481 eV for case (b). 


In order to investigate birefringence effects in the active 
region of a conventional laser, we consider uniaxial strain 
by extending the lattice constant in x-direction. For sim¬ 
plicity, we assume the barrier to have the same lattice 
constant as GaAs, 5.6533 A, in y-direction. Therefore, 
both barrier and well regions will have the same exten¬ 
sion in x-direction. For Qx = 5.6544 A we have the cor¬ 
responding element of the strain tensor Sxx ^ 0.019%, 
while ttx = 5.6566 A gives £xx ~ 0.058%. 

The effect of uniaxial strain in the band structure is 
presented in Fig. dHa) and (b) for £xx ^ 0.019% and 
£xx ^ 0.058%, respectively. The labeling and ordering of 
subbands follows the same as that in Fig. [Sl]a). Just this 
slight anisotropy in the x- and y-lattice constants creates 
a difference in subbands for the [ 100 ] and [ 010 ] directions. 
In the inset we show the region around the anti-crossing 
of HHl and LHl subbands, where the difference is more 
visible. 

In addition to the differences induced in the band 
structure, the uniaxial strain also induces a change in 
the dipole selection rules between and S'^ light po¬ 


larizations, which can be seen in the gain spectra we 
present in Fig. fT^ a) and fl^ b) for Sxx ~ 0.019% and 
£xx ~ 0.058%, respectively. Reflecting the features of the 
band structure, we notice for the emission region of the 
gain spectra that the largest difference between and 
gy is around the HHl and LHl energy regions (between 
Cl and C 3 cavity positions, approximately). In the ab¬ 
sorption regime (negative gain) we notice g^ < g^, while 
in the emission regime (positive gain) we have g^ > g^. 
This feature is more visible in Fig, fl^ b). 


C2 C3 



Eph - Eg (meV) 

Figure 12: (Color online) Uniaxial strain modification of gain 
spectra for strain (a) £x;x ~ 0.019% and (b) e^x ~ 0.058%. 
The anisotropy in the lattice constants for x- and y-directions 
modifies the output light polarization of the laser. Since there 
are no spin-polarized carriers in the system, ~ 9 - 


To calculate the birefringence coefficient in the active 
region, we used the definition of Ref. [ 6 ^ given by 

7p(w) = ^fer(w) , (6) 

ZTleilg 

where uj is the frequency of the longitudinal mode in the 
cavity, rig is the effective index of refraction of the cavity, 
and rig is the group refractive index. For simplicity, we 
assume rig = rig. The real part of the dielectric func¬ 
tion can be obtained from the imaginary part using the 
Kramers-Kronig relations4^ 

We present the birefringence coefficient in Fig. fT^ a) 
and 113Kb) for £xx ^ 0.019% and £xx 0.058%, re¬ 
spectively. We notice that this strain in the active re¬ 
gion, responsible for modest changes in the gain spec¬ 
tra, produces birefringence values of the order of 10 ^^“^^ 
Hz which may be exploited to generate fast polariza¬ 
tion oscillations. Furthermore, when increasing the strain 
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amount by ^ 0.04% from case (a) to case (b), the value 
of 7 p increases approximately threefoldi^ We also in¬ 
cluded in our calculations spin-polarized electrons and 
notice that they have only a small influence in the bire¬ 
fringence coefficient. Although they change \g^\ and \gy\ 
slightly , the asymmetry is not affected at all for small 
spin polarizations of 10 - 20 %, which are relevant values in 
real devices. 



Eph - Eg (meV) 

Figure 13: (Color online) Birefringence coefficient as a func¬ 
tion of photon energy considering (a) ~ 0.019% and 

(b) Exx ~ 0.058%. Just an increase of 0.0022 A in Ux in¬ 
creases 7 p by approximately 3 times. The two peaks, around 
Eph — Eg ~ 0 meV and Eph — Eg ~ 150 meV are related to 
transitions from CBl and CB2. Transitions related to CB2 
are in the absorption regime, not visible in Fig. 1121 


Investigating the effect of different cavity designs, we 
present the values of 7 p in Figs. fT^ a) and fT^ b) for 
Exx ^ 0.019% and Sxx ~ 0.058%, respectively. We chose 
the same photon energies as for the case without birefrin¬ 
gence assuming that the different values for the strain- 
induced birefringence in the active region will not signif¬ 
icantly affect the cavity resonance for reasons of simplic¬ 
ity. For the two different strain types the behavior of ■jp 
is very similar for the same resonance energy. Comparing 
different cavity designs we observe that for Ci, the value 
of 7 p strongly decreases and also changes sign with the 
carrier density, n. In contrast, for C 2 and C 3 , ■jp is always 
positive. After a slow increase with n, 7 ^ becomes flat, 
and nearly independent of the carrier density. 

For consistency, we have also calculated the DBR con¬ 
tributions using the approach given by Mulet and Ballej^ 
For large anisotropies in the DBR, the birefringence co¬ 
efficient is on the order of 10^° Hz, consistent with the 
measurements given by van Exter et ali^ Therefore, for 



Figure 14: (Color online) Birefringence coefficient as function 
of the carrier density for (a) Exx 0.019% and (b) Exx ~ 
0.058%. For different cavity designs the behavior of 7 ^ can 
be completely different. The carrier density values where 7 p 
changes sign in cavity ci and the flat region in cavities C 2 and 
C 3 are already in the lasing regime. 


the investigated strain conditions, the main contribution 
to 7 p comes from the active region and it is a very versa¬ 
tile parameter that can be fine-tuned using both carrier 
density and cavity designs, possibly even changing its 
sign and reaching carrier density-independent regions. 

VII. ULTRA HIGH-FREQUENCY OPERATION 

Lasers could provide the next generation of par¬ 
allel optical interconnects and optical information 
processingThe growth in communication^ 
and massive data centers^ will pose further limitations 
on interconnects!^ Conventional metallic interconnects 
used in multicore microprocessors are increasingly recog¬ 
nized as the bottleneck in maintaining Moore’s law scal¬ 
ing and the main source of power dissipation.—!^ Optical 
interconnects can effectively address the related limita¬ 
tions, such as the electromagnetic crosstalk and signal 
distortion, while providing a much larger bandwidth— 
VCSELs are considered particularly suitable for short- 
haul communication and on-chip interconnects— How¬ 
ever, to fully utilize their potential, it would be impor¬ 
tant to explore the paths for their high-frequency opera¬ 
tion and achieve a higher modulation bandwidth, limited 
for conventional lasers to about ^ 50 GHz.— 

How can we understand the frequency limitation of a 
laser? Why would a higher frequency modulation lead 
to a decrease in a signal-to-noise ratio and limit the ef¬ 
fective bandwidth? An accurate analogy is provided by 
a driven and damped harmonic oscillator. The laser re¬ 
sponse, just like the harmonic oscillator, is unable to fol¬ 
low a high enough modulation frequency. A Lorenzian- 
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like frequency-dependent displacement of a harmonic os¬ 
cillator closely matches a modulation response of a laser, 
decreasing as l/w, above the corresponding resonance 
frequency, known as the relaxation oscillation frequency, 
Wfl, representing a natural oscillation between the car¬ 
riers and photons and often used to estimate the band¬ 
width of a laser 

To realize a high-speed operation in conventional lasers 
requires a careful design and optimization of many pa¬ 
rameters. Attaining a high w/j is closely related to opti¬ 
mizing the gain which increases with but decreases 
with photon density S', known as the gain compression^^ 
which would be desirable to minimize. For a small-signal 
modulation S{t) = Sq + SS{t), above the threshold,— 

Vg{dg/dn)So/TpYi, (7) 

where Vg is the group velocity of the relevant mode, 
dg/dn is the differential gain at the threshold, and Tph 
is the photon lifetime. While ujr increases with So, a 
larger Sq, through gain compression, is detrimental by 
diminishing the differential gain. There are additional 
factors, beyond Eq. 0 , required for a high w/j, such as 
minimizing the transport time for carriers to reach the ac¬ 
tive region, achieving a high carrier escape rate into the 
QW barriers, and minimizing extrinsic parasitic effects 
between the intrinsic laser and the driving circuit 

Introducing spin-polarized carriers offers additional 
possibilities to enhance w/j, corresponding to the modu¬ 
lation of the emitted S, beyond the frequencies attainable 
in conventional lasers. In the regime of small-signal mod¬ 
ulation, both uiR and the bandwidth have been shown to 
increase with an increase of the spin-polarization of the 
injected carriers, associated with the threshold 

reduction [thus for a given injection So is larger than 
in Eq. ©]• Similar trends are predicted in the large- 
signal modulation, but the corresponding increase of lor 
(as compared to the conventional lasers) can exceed what 
would be expected based only on the threshold reduction 
due to Pj ^ 0— 

Another approach to achieve a higher ujr is to use the 
polarization dynamics, instead of the intensity dynamics 
of the emitted light. The coupling between spin-polarized 
carriers and the light polarization in birefringent mi¬ 
crocavities corresponds to different resonant mechanisms 
than those that govern the light intensity and thus to 
potentially higher lur. Early experiments on polariza¬ 
tion dynamics in VCSELs of Oestreich and collaborators 
have demonstrated spin-carrier dynamics of 120 GHz.— 
However, their (Ga,In)As QW spin lasers operated at 10 
K and required a large magnetic field for fast spin pre¬ 
cession. 

Gould we attain similar ultrahigh frequencies at room 
temperature without an applied magnetic field? Our 
findings from Sec. |Vl] suggest that indeed such an op¬ 
eration could be realized by a careful design of bire¬ 
fringent cavity properties providing frequency splitting 
of the two orthogonal linearly-polarized lasing modes. 


While in conventional VC SEE only one linearly-polarized 
mode is emitted, injecting spin-polarized carriers leads 
to the circularly-polarized emission and thus the opera¬ 
tion of both linearly-polarized modes at the same time. 
The beating between the two frequency-split linearly- 
polarized modes creates polarization oscillations with fre¬ 
quency determined by the birefringence rate, 7 p/ 7 r.— 
Strain-induced values of jp in the active region shown 
in Figs. [13] and [T3| are sufficiently high to exceed the 
highest available frequency operation of conventional VC¬ 
SELs. A strong spectral dependence of 7 p, including a 
possible sign change, requires a careful analysis of the 
detuning behavior, but it also provides important op¬ 
portunities for desirable operation of spin lasers. For 
example, a large 7 p can be achieved with a very weak 
dependence on the carrier density. The feasibility of 
a high-birefringence rate is further corroborated by the 
experiments using mechanical strain attaining Jp/tt ^ 
80 GHz,— while theoretical calculations suggest even 
7 p/ 7 r ^ 400 GHz with asymmetric photonic crystals— 


VIII. CONCLUSIONS 

Our microscopic model of optical gain is based on a 
similar framework previously employed for conventional 
lasersdLSidi to simply elucidate how introducing spin- 
imbalance could enable their improved dynamical opera¬ 
tion. In contrast to the common understanding that the 
birefringence is detrimental for lasers, we focus on the 
regime of a large strain-induced birefringence to over¬ 
come frequency limitations in conventional lasers. 

With a goal to maximize the birefringence-dominated 
bandwidth in a experimentally realized spin laser, 
we can use the guidance from the analysis of both 
high-speed conventional lasers and the steady-state 
operation of spin lasers to explore potential limit¬ 
ing factors. Future calculations should also examine 
the influence of a spin-dependent gain compression. 
Coulomb interactions an active region with multi¬ 

ple QWs^, spin relaxation^SiS^iZ^ and a careful analysis 
of the optimal cavity position that would combine high 
(differential) gain, high-gain asymmetry,and high 7 ^. 

While currently the most promising path to demon¬ 
strate our predictions for ultrahigh frequency operation 
is provided by optically injected spin-polarized carriers 
to the existing VCSELs, there are encouraging devel¬ 
opments for electrically injected spin-polarized carriers. 
A challenge is to overcome a relatively large separation 
between a ferromagnetic spin injector and an active re¬ 
gion (> /im) implying that at 300 K recombining car¬ 
riers would have only a negligible spin polarization— 
However, room temperature electrical injection of spin- 
polarized carriers has already been realized through spin¬ 
filtering by integrating nanomagnets with the active re¬ 
gion of a VCSEL— Additional efforts focus on vertical 
external cavity surface emitting lasers (VECSELs)^d^ 
which could enable depositing a thin-fihn ferromagnet 
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to be deposited just 100-200 nm away from the active 
region, sufficiently close to attain a considerable spin po¬ 
larization of carriers in the active region at room temper¬ 
ature. 

An independent progress in spintronics to store and 
sense information using magnets with a perpendicular 
anisotropy^ and to attain fast magnetization reversal^ 
could also be directly beneficial for spin lasers. Electri¬ 
cal spin injection usually relies on magnetic thin films 
with in-plane anisotropy requiring a large applied mag¬ 
netic field to achieve an out-of-plane magnetization and 
the projection of injected spin compatible with the car¬ 
rier recombination of circularly polarized light in a VC- 
SEL geometry (along the z-axis, see Fig. |4]). How¬ 
ever, a perpendicular anisotropy could provide an elegant 
spin injection in remanence^^— avoiding the technologi¬ 
cally undesirable applied magnetic field. The progress in 
fast magnetization reversal could stimulate implement¬ 
ing all-electrical schemes for spin modulation in lasers 
that were shown to yield an enhanced bandwidth in 
lasers 

Note added in proof. After this work was completed 
and submitted, our predictions for high-frequency bire¬ 
fringence were experimentally demonstrated in similar 
GaAs/AlGaAs quantum well spin VCSELs revealing val¬ 
ues of ^ 250 GHa22. 


Acknowledgements 

We thank M. R. Hofmann and R. Michalzik for valu¬ 
able discussions about the feasibility of the proposed spin 
lasers and the state-of-the art active regions in conven¬ 
tional lasers. We thank B. Scharf for carefully read¬ 
ing this manuscript. This work has been supported by 
GNPq (grant No. 246549/2012-2), FAPESP (grants No. 
2011/19333-4, No. 2012/05618-0 and No. 2013/23393- 
8 ), NSF EGGS-1508873, NSF ECCS-1102092, U.S. ONR 
N000141310754, NSF DMR-1124601, and the German 
Research Foundation (DFG) grant ’’Ultrafast Spin Lasers 
for Modulation Frequencies in the 100 GHz Range” GE 
1231/2-1. 


Appendix A 


lasers as well as to elucidate a wealth of other 

phenomena, such as the spin Hall effect, topological in¬ 
sulators, and Zitterbewegung.— Our own implementa¬ 
tion of the k-p method in this work has been previously 
tested in calculating the luminescence spectra in (5-doped 
GaAsfSi confirming experimental and theoretical elec¬ 
tronic structure for GaAs QWs,— and (Al,Ga)N/GaN 
superlattices^ identifying fully spin-polarized semicon¬ 
ductor heterostructures, based on (Zn,Go)Ori and ex¬ 
ploring polytypic systems consisting of zinc-blende and 
wurtzite crystal phases in the same nanostructure.— > 2 ^ 

Before considering confined systems, it is important to 
investigate the corresponding bulk crystal structure and 
construct the functional form of the Hamiltonian. For 
zinc-blende crystals, the bulk basis set that describes the 
lower conduction and top valence bands igSSnSiviflS 


|CB^) = |5t) 

|CB^) = |5|) 

|HH^) = \{X + iY)f)/y/2 
|LH^) = i\{X + iY)i-2Zf)ly/Q 
|LH^) = \{X-iY)\+2Zi)/^f& 

|HH4) = i\{X-iY) i)/V2 
|SO^) = \{X + ^Y)l+Z\) 

|SO^) = i\-{X-iY)\+Z i) /y/i, (A-1) 

where, compared to Fig. HKa), we also introduce the 
spin-orbit spin-split-off subbands |SO). Here [S') and 
I A), |y), \ Z) are the basis states for irreducible repre¬ 
sentations Pi ^ and Pis ~ x^y,z, having 

an orbital angular momentum I = 0 and Z = 1 , respec¬ 
tively. The single arrows (t,/-) represent the projection 
of spin angular momentum s = 1/2 on the -I- z-axis while 
the double arrows (jf, IJ.) represent the projection of total 
angular momentum on the -I- z-axis. Rewriting the basis 
set (lA-ll) in terms of the total angular momentum j and 
its projection rrij, \j,mj), we have 

|GB^(^)) = |l/2, 1/2 (-1/2)) 

|HH^(^)) = |3/2, 3/2 (-3/2)) 

|LH^(^)) = |3/2, 1/2 (-1/2)) 

|SO^(^)) = |l/2, 1/2 (-1/2)) . (A-2) 


The versatility of the k-p method has been success¬ 
fully used to obtain the gain spectra in conventional 


In the basis set of Fa. (lA-ll). the k-v term in Fa. is 
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(A-3) 


with elements 

Q = -kx (71 + 72) kx - ky (71 + 72) ky - k^ (71 - 272) k^ 
T = -kx (71 - 72) kx - ky (71 - 72) ky - k^ (71 + 272) k^ 

S = iVi[{kx^skz + kzlzkx)-i{ky^:ikz+k:^^:iky)\ 

R = [{kx^2kx - kyj2ky) - i {kxjsky + kyjskx)] 

U = Eg kx-^kx kyAky kz-Akz 

P± = {I/ 2 V 2 ) [P {kx ± iky) + {kx ± iky) P] 

P, = (1/2) {Pk, + k,P) , (A-4) 

where 71, 72,73, and A, given in units of h'^/ 2 mo, are the 
effective mass parameters of the valence and conduction 
bands, respectively, explicitly given below. The gap is 
Eg, the spin-orbit splitting at the T point is A50, and 
P is the Kane parameter of the interband interaction, 
defined as 


P =-i—{a\pi\S) , (A-5) 

Too 

with a = X,Y, Z and I = x,y, z. 

The formulation of a bulk k-p model can vary signifi¬ 
cantly in its complexity, the choice of the specific system, 
and the number of bands included. In the description of 
zinc-blende structures, usually either 6x6 or 8x8 models 
are employedi^S In the first case, the information of the 
valence and conduction band is decoupled, while in the 
second case their coupling is explicitly included. Their 
effective mass parameters are connected by 

71 = 71 - Ep/3Eg 

72 = 'y2-Ep/QEg 

73 = 'js-Ep/QEg 

^ ^ _ f Eg + 1^50 \ Ep 

\ Eg + Aso ) Eg 

Ep = 2moPV^^ (A-6) 

where 71,2,3 are used in the 8x8 model and 71,2,3 in the 
6x6 model, which can also be related to the tight-binding 


parameters.— To recover the 6x6 model from the 8x8 
model, we set P = 0 in Eqs. (IA-31) . (IA-41) and (IA-6|) . 

The strain term, Hgt, takes a form similar to Eq. (IA-31) 
but without the Eg, Aso and P parameters. The matrix 
elements can be written as 

Qst — 

Tst = 

Sst = 

Rst = 

Ust = 

with tty, b, and d representing the deformation potentials 
for the valence band and Qc for the conduction band. The 
strain tensor components are given by {i,j = x,y,z). 

In order to treat a QW system, which now lacks trans¬ 
lational symmetry along the growth direction, we can 
replace the exponential part of the Bloch’s theorem by a 
generic function. This procedure is called the envelope 
function approximation^ and it leads to the dependence 
along the growth direction of the k-p and strain parame¬ 
ters in Hamiltonian terms Hkp{z) and E[st{z). Also, the 
band-offset at the interface of different materials is taken 
into account in the term E[o{z) 

Ho{z) = diag [dv{z), ■■■ , 6v{z), Sc{z), 5c{z)] , (A-8) 

where 5v{c) {z) describes the energy change in the valence 
(conduction) band. 

Under the envelope function approximation, the QW 
Hamiltonian from Eq. © is now described by a system 
of 8 coupled differential equations that does not gen¬ 
erally have analytical solutions. We solve these equa¬ 
tions numerically using the plane-wave expansion for the 
z-dependent parameters and envelope functions. De¬ 
tails of the envelope function approximation and plane 
wave expan sion for QW systems can be found in refer¬ 
ences 


i^^xx ^yy E ^zz) 2 ^^xx ^yy ^^zz') 

O'v i^^xx ^yy E ^zz') {^xx ^yy ^^zz') 

d {f^yz “ 1 “ ^'^ 3 : 2 :) 

_ \\ide 

2 \^xx ^yyJ ' f'^^xy 
i^^xx ^yy E ^zz) j 
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Appendix B 


and therefore 


The interband dipole transition amplitude that ap¬ 
pears in Eq. m is given by 


a ■ p = 


Px ± ipy 


(B-3) 


P^vk = {c,k\a-p\v,k) , 


and for the light polarization we have 
a= ^{x±iy) , 


In the simplified QW of Fig. 01 we are showing the 
selection rules for = 0 and assuming the conduction 
band as |c, 0) = |CB (H)), and valence band as, |u, 0) = 
|HH fl" (J|)) or |u, 0) = |LH (IJ.)). Calculating the matrix 
(B-2) elements between these states, we obtain 


(CB^|p±|HH^) = (5t 

which is non-zero only for p_ 

(CB^|p±|HH^) = 
which is non-zero only for p_|_, 
(CB^|p±|LH^) = /^t 


Px ± iPy 


V2 

Px ± iPy 


72 


(X + iY) t 


Px ± iPy 


72 

which is non-zero only for p-|_, and 

Px ± iPy 


72 


1 




— [{X-xY)t+2Zi] 


(CB^|p±|LH^) = (^,5; 

which is non-zero only for p_. 


72 


76 


[(A + ^Y) i -2Z t] 


^{St\Px\Xt)Tl{St\Py\Yt) , (B-4) 

'-{Si\p,\Xi)±^-{Si\py\Yi) , (B-5) 

:^{St\Px\Xt)±^{St \Py\Y t) (B-6) 
:^{Si\p,\Xi}T^{Si \py\Y i) (B-7) 


In addition to Eqs. (IB-4ll - (|B-7(l . we can conclude (CB ff IpilLHlI) = (CB |p±|LHJJ.) = 0, independent 

that (CB fl" |p± I HH H) = (CB IJ. |p±| HH fl-) = 0 and of the light polarization. 
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